The reconstruction of Ni and Rh (001) surfaces 
upon Carbon, Nitrogen, or Oxygen adsorption 



Dario Alfe,B'@ Stefano de Gironcoli,B and Stefano BaroniB 
SISSA - Scuola Internazionale Superiore di Studi Avanzati and 
INFM - Istituto Nazionale di Fisica della Materia 
via Beirut 2-4, 1-34014 Trieste, Italy 
(February 1, 2008) 

Abstract 

Nickel and Rhodium (001) surfaces display a similar — as from STM images — 
clock reconstruction when half a monolayer of C/Ni, N/Ni or O/Rh is ad- 
sorbed; no reconstruction is observed instead for O/Ni. Adsorbate atoms 
sit at the center of the black squares of a chess-board, c(2 x 2), pattern and 
two different reconstructions are actually compatible with the observed STM 
images — showing a (2 x 2)p4g pattern — according to whether a rotation of 
the black or white squares occurs. We report on a first-principles study of the 
structure of X/Ni(001) and X/Rh(001) surfaces (X=C,N,0) at half a mono- 
layer coverage, performed using density-functional theory. Our findings are in 
agreement with all available experimental information and shed new light on 
the mechanisms responsible for the reconstructions. We show that the same 
substrate may display different reconstructions — or no reconstruction — upon 
adsorption of different atomic species, depending on the relative importance 

of the chemical and steric factors which determine the reconstruction. 
PACS numbers: 68.35.Bs 82.65.My 82.65.Jv 
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I. INTRODUCTION 



Transition metal surfaces are attracting a wide scientific and technological interest, es- 
pecially because of their capability of reducing the activation barrier of many important 
chemical reactions. One particularly relevant example is the reaction 2CO + 2NO — > 2CO2 
+ N 2 , which eliminates the two poisonous CO and NO gases from the exhaust gas of com- 
bustion engines. This reaction is catalysed by many transition metal surfaces, among which 
rhodium and platinum have been demonstrated to be the most efficient. 

Adsorption of atomic or molecular species on surfaces modifies their electronic and struc- 
tural properties, thus affecting their catalytic properties. Recently, the (001) surfaces of 
rhodium and nickel have been studied with respect to their peculiar reconstructions upon 
oxygen, carbon and nitrogen adsorption [fl| — 110 1 . Oxygen adsorption on Rh (001) is known 
to be dissociative and to saturate at half a monolayer, independently of the adsorption 
temperature. At this coverage a (2 x 2) reconstruction has been observed by LEED || 
and confirmed by STM ||. The oxygen atoms sit in the troughs formed by four first-layer 
rhodium atoms and fill these sites in a c(2 x 2) geometry. This structure may be seen as 
a chess-board whose black squares are occupied by oxygen atoms, while the white ones are 
empty. Within this picture, the reconstruction observed in Ref. [[J has been described as 
a rotation of the black squares, resulting in a (2 x 2)p4g symmetry (see Fig. |T]a) . This 
distortion preserves the shape of the black squares, while the white ones become rhomboid. 
A similar behavior is observed for nitrogen and carbon adsorbed on the (001) surface of 
nickel, where the rotation angle of the squares is much larger and the clock reconstruction 
more evident |5|-|10"1. 

In our previous work on Rh(001)/O @] we pointed out that a different substrate re- 
construction is actually compatible with LEED and STM data. Using only these two ex- 
perimental techniques it is not possible to distinguish between the reconstruction described 
above and the one which results instead from the rotation of the white squares (Fig. |l]b). 
Our calculations actually indicated that it is the white squares which rotate rather than 
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(a) (b) 

FIG. 1. Two possible clock reconstructions of the O/Rh(001) surface, resulting in a (2 x 2)p4g 
structure. In the black reconstruction (a) the squares with an O atom in the middle rotate, while 
the empty ones distort to rhombi; in the white reconstruction (b), the opposite occurs. 

the black ones. Moreover, we predicted a reconstruction pattern with a slightly different 
symmetry from what appears in the STM pictures. We found that the ad-atoms get off 
the center of the distorted squares, thus resulting in an asymmetric clock reconstruction 
(Fig. ||). The oxygen atoms can occupy two equivalent low-symmetry sites separated by 
a low energy barrier, of the order of the room thermal energy. This would give rise to an 
order-disorder transition, and we suggested that what is actually observed in the STM pic- 
tures is the average position of the oxygen atoms jumping back and forth between the two 
equilibrium positions, in the high-temperature disordered phase. 

The reconstruction of the (001) surface of nickel upon C and N adsorption is experimen- 
tally much better characterized since also LEED-IV || and SEXAFS || data exist. In this 
case it is widely accepted that the reconstruction is of the black type, i.e. it is the filled 
squares which rotate (Fig. |I|a). No surface reconstruction is induced by oxygen adsorption 
on Ni(001) [jnj, while as far as we know no experimental data exist for carbon or nitrogen 
adsorbed on the Rh(001) surface. 

In this paper we present an ab-initio study of the adsorption of carbon, oxygen, and 
nitrogen on the (001) surfaces of rhodium and nickel. On rhodium, we find that neither 
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FIG. 2. Equilibrium structure of O/Rh(001) as obtained by our simulated-annealing procedure. 
The thick line indicates the unit cell. The thin line indicates the rhombus and its shorter diagonal. 
The oxygen atoms are alternatively shifted orthogonally with respect to the shorter diagonal. The 
atoms of the first surface layer depicted with a brighter tone lean ~ 0.08 A outwards with respect 
to the others. 

nitrogen, nor carbon, induce any reconstruction of the surface. On nickel, in agreement with 
the experiments, we find that both carbon and nitrogen induce a clock reconstruction, while 
oxygen adsorbs with no induced reconstruction. We discuss the interplay of the chemical 
and the steric effects which is at the root of these behaviors. 

The paper is organized as follows. In section we describe our theoretical method, 
including some tests on the properties of the rhodium and nickel bulk metals, and on those 



of the CO, 2 , NO and N 2 molecules. In section |T| we present our results. Section [TV 
contains the discussion and our conclusions. 

II. METHOD 

Our calculations are based on density functional theory within the local-density approx- 
imation (LDA) |12[], using Ceperley-Alder exchange-correlation (XC) energies The 
one-particle Kohn-Sham equations are solved self-consistently using plane- wave basis sets in 
a ultra-soft (US) pseudo-potential scheme [inj. The rhodium, oxygen, and carbon pseudo- 
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potentials are the same as in Refs. [|T^,|T^]. For nitrogen and nickel we have constructed new 
pseudo-potentials. In the case of Ni, we started from the 3<i 8 4s 2 4p° reference configuration 
treating the d channel in the US scheme, while the s and the p channels were assumed to 
be norm conserving. We chose the p component of the pseudo-potential as local part (this 
choice avoids the appearance of ghost states). For N we started from the 2s 2 2p 3 reference 
configuration treating both the s and the p channels in the US scheme. Plane waves up 
to maximum kinetic energy of 30 Ry are included in the basis set. Brillouin-zone integra- 
tions have been performed using the Gaussian-smearing fl7| special-point fl8| technique. 



For the calculations in the bulk, we have used a smearing function of order 1 with a width 
a = 0.03 Ry and a set of special points equivalent to 10 k-points in the irreducible wedge 
of the Brillouin zone (IBZ). The convergence of our results with respect to the size of the 
basis and special-point sets has been carefully checked. 

The isolated surfaces are modeled by a periodically repeated super-cell. We have used 
the same super-cell for both the clean and the covered surfaces. For the clean surfaces we 
have used 5 atomic layers plus a vacuum region corresponding to 6 layers. For the covered 
surfaces, the 5 Rh/Ni layers are completed by one layer of C, O, or N atoms on each side of 
the slab; in this case the vacuum region is correspondingly reduced to ~ 5 atomic layers. We 
have used a Gaussian-smearing function of width o = 0.03 Ry and a (12, 12, 2) Monkhorst- 
Pack mesh |18| resulting in 21 special k-points in the lxl surface IBZ. Convergence tests 



performed with a value of a twice as small and a correspondingly finer mesh of special points 
resulted in no significant changes in total energies and equilibrium geometries. 

A. Case tests on bulk metals and molecules 

In order to test the quality of our pseudo-potentials, we have done some tests on rhodium 
and nickel bulks, as well as on CO, NO, O2, and N2 molecules. To calculate the properties of 
the molecules the latter have been put in a large cubic cell, with side L = 10 a.u, and periodic 
boundary conditions. We have checked that our calculated molecular properties were well 
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converged with respect to the size of the cell. The results of these tests are summarized in 



Tab. | and Tab. |TJ. 

TABLE I. Theoretical structural parameters of FCC rhodium and nickel obtained with a fit to 
Murnaghan equation of state of the calculated energies as a function of the volume. Experimental 



data [19 are reported for comparison. 



Lattice constant Bulk modulus 

(A) (Mbar) 



Rh Present work 3.81 3.17 

Experiment 3.80 2.69 

Ni Present work 3.44 2.4 

Experiment 3.52 1.88 



TABLE II. Theoretical equilibrium distances and fundamental vibrational frequencies. Ex- 
perimental data pH] are reported for comparison. 







Equilibrium distance 

(A) 


Vibrational frequency 
(cm -1 ) 


CO 


Present work 


1.14 


2200 




Experiment 


1.13 


2170 


NO 


Present work 


1.16 


2000 




Experiment 


1.15 


1904 




Present work 


1.23 


1650 




Experiment 


1.21 


1580 


N 2 


Present work 


1.11 


2450 




Experiment 


1.09 


2360 



Nickel and rhodium have the FCC crystal structure. To calculate the equilibrium lattice 
constant ao and the bulk modulus B , we have fitted the calculated energies as a function 
of the unit cell volume to a Murnaghan's equation of state. Nickel is a magnetic metal 
with a magnetic moment of 0.59 /is/atom [ffi}| . An explicit account of the spin polarization 
within the local spin density approximation did not result in any meaningful changes of 
the calculated structural properties. The magnetism is expected to be more important at 
the surface, because the surface density of states (SDOS) is narrower than the bulk one. 
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FIG. 3. d-projected density of states of nickel for the bulk, the clean surface, and the car- 
bon-covered surface. 

However, at least for the carbon- and the nitrogen-covered surfaces, we have found that the 
Ni-<i-SDOS is shifted towards lower energies with respect to the bulk d partial density of 
states (see Fig. ^) thus resulting in a lower density of states at the Fermi level. We argue 
that magnetic effects should not be very important and they have been neglected altogether 
in the present investigation. 



III. RESULTS 

In order to find the stable surface structure of the C-, 0-, and N-covered Rh(001) and 
Ni(001) surfaces, we have performed simulated annealing runs of the surfaces at half a 
monolayer of coverage, using Born-Oppenheimer ab-initio molecular dynamics (AIMD). In 
our implementation of AIMD the exact electronic ground state (within a self-consistent 
threshold) is calculated at each time step using a band- by-band conjugate gradient algorithm 
|22|| and Broyden-like [23] mixing of the potentials. The forces are calculated fully quantum 



mechanically, and the ionic equation of motion is integrated using the Verlet algorithm p3 



TABLE III. Summary of the structural data for the six systems investigated in this work. For 
the reconstruction type see Fig. [l|. doi is the distance between the adatoms and the first metal layer; 
S is the amplitude of the in-plane displacement of the first-layer metal atoms upon reconstruction; 
AE is the energy difference between the symmetric c(2 x 2) structure and the reconstructed one. 
Units are A for lengths and eV/atom for energies 





Reconstruction 
Expt Theory 


5 

Expt 


Theory 


doi 

Expt 


Theory 


AE 
Theory 


C:Ni 


black 


black 


0.55±0.20 a 


0.46 


0.1±0.1 a 


0.17 


0.20 


N:Ni 


black 


black 


0.55 ± 0.20 a 


0.41 


0.1±0.1 a 


0.10 


0.08 


0:Ni 


none 


none 






0.77±0.04 b 


0.73 




C:Rh 




none 








0.63 




N:Rh 




none 








0.71 




0:Rh 


white 


white 


0.20 ± 0.07 cd 


0.21 


1.05±0.05 c 0.6±0.1 d 


0.98 - 1.06° 


0.03 



a From Ref. [1C] and references quoted therein 
b From Ref. ITT | 



c From Ref. (25 1 
d See Ref. @ 

e The two values correspond to the two Rh rows made inequivalent by the asymmetric clock recon- 
struction (see text and Fig. |^). 

We have used a time step of 50 a.u (~ 1.2 fs). After an equilibration period of about 0.5 
ps at a temperature of ~ 800 K we have slowly cooled down the systems. Finally, when 
the ionic minimum structure is reached, an accurate structure optimization is performed, by 
allowing all the atoms in the slab to relax until the force acting on each of them is smaller 
than 0.5 x 10 _3 Ry/ao. Our results are summarized in Table JTJ and discussed in detail in 
the following sections. 



A. Rhodium 

Results for the O/Rh(001) surface have been reported and extensively discussed in our 
previous work ||]. For this system we have found that oxygen induces an asymmetric clock 
reconstruction of white type (Fig. |2|), where the oxygen- filled sites are deformed into rhombi, 
and the ad-atoms get off the center of the rhombi. The asymmetric position of the oxygen 
ad-atoms results in a buckling reconstruction of the first surface layer, and the rhodium rows 
are alternatively shifted up and down by ~ 0.08 A. Oxygen atoms tend to shorten the bonds 
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with the neighboring rhodium atoms, and this results in a deformation of the oxygen site 
which becomes a rhombus. We argued that the driving mechanism for this reconstruction 
is the re-bonding of the oxygen with the first layer rhodium surface. This can be better 
understood by the inspection of the local density of states (LDOS). 

The LDOS is the projection of the density of states onto localized atomic orbitals. In Fig. 
^a we display the Kh(4d xz ) LDOS of one rhodium surface atom resulting from a calculation 
with the clean surface; this is degenerate with the Kh(4d yz ) orbital because of symmetry. 
The energy is referred to the Fermi level. In the panel just below (b) we display the same 
two orbitals after half a monolayer of oxygen is deposited, without allowing the substrate to 
reconstruct. The Kh(4d xz ) orbital points towards the square site left empty, and its LDOS 
is very similar to what it would be for the clean surface. The Kh(4d yz ) orbital is directed 
towards the oxygen adsorption site, and the shape of its LDOS changes completely upon 
adsorption: the most important feature is the presence of two peaks, one at pa —6 eV and the 
second above the Fermi energy, at ~ 1 eV. In the (c) panel we display the LDOS projected 
onto 0{2p x ) orbitals, which are degenerate with the 0(2p y ) orbitals because of the symmetry 
of the oxygen site. This LDOS is characterized by a bonding and an anti-bonding peak, 
at energies of ~ —6 eV and pa +1 eV respectively, which are in the same positions of the 
peaks observed in the Rh(4d yz ) LDOS. We interpret this fact as the evidence of the presence 
of a strong covalent bond between the oxygen atoms and the rhodium surface atoms. The 
different weight between the bonding and the anti-bonding peaks also indicates that there 
is a charge transfer from the substrate to the oxygen, so that the bond is partially ionic. 
This fact is also supported by work function calculations, which show an increase of ~ 0.6 
eV when the surface is covered by oxygen. This imply that an excess of negative charge 
is present on the surface. In the last panel of the figure we display the LDOS projected 
onto the 0(2p x ) and 0(2p y ) orbitals for the reconstructed surface. The reconstruction lifts 
the xy symmetry, and the 0(2p x ) orbital — which corresponds to the shorter O-Rh bond — 
becomes more occupied than the 0(2p y ) orbital. We conclude that the mechanism of the 
reconstruction is due to the re-bonding of the oxygen with the surface atoms, which tend 
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c(2x2), Oxygen p x = p y 




-8 -7 -6 -5 -4 -3-2-1012 




-8 -7 -6 -5 -4 -3-2-1012 



FIG. 4. Local density of states projected onto surface atomic orbitals. (a): Rh d xz and d yz 
orbitals of the clean surface, (b): Same as (a), but after the deposition of half a mono-layer of 
oxygen, without allowing the substrate to reconstruct, (c) Oxygen p x and p y orbitals. (d) Same 
as (c), but after the clock reconstruction has taken place. 
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to shorten the O-Rh bond. This shortening could also be realized by a penetration of the 
oxygen deeper into the site, but this does not happen. The reason is maybe the excess of 
negative charge on the oxygen atoms, which prevents their penetration into the electronic 
sea of the metallic surface. 

While the present work was being completed, we learnt that our theoretical predictions 
have been recently confirmed by LEED I-V experiments [P5|| . In particular, the analysis 
of the LEED I-V spectra indicated that the reconstruction of the Rh surface upon oxygen 
adsorption is indeed of the white type. Furthermore, the stable adsorption site of the oxygen 
adatoms was found to be asymmetric with respect to the center of the rhombus, with a 
displacement 5o = 0.29 ± 0.15 A, to be compared with our theoretical prediction, <5 G = 
0.35 A. These data are in partial agreement with the recently published results of an 
experiment based on a combination of LEED with low-energy alkali ion scattering and 
recoil spectroscopy [ 26| . In Ref . [26] no evidence was found of the asymmetric white clock 



reconstruction predicted by us and confirmed in Ref. f25(| . Furthermore, the height of the 
O adsorption site off the metal surface was found to be smaller than our predictions which, 
however, are again in agreement with the findings of Ref. EH . 



At variance with oxygen, nitrogen does not induce any reconstruction upon adsorption 
on Rh(001). The equilibrium distance of the adatom from the metal surface is rather large, 
as it is the case for oxygen (see Table For carbon we have found a metastable state with 
a surface reconstruction of the black type (Fig. |l]a) and a rather small distance between the 
ad-atoms and the metal surfaces (see the discussion of the black reconstructions typical of 
Ni in the next subsession) . However, the energy of this local minimum is higher than that 
of the unreconstructed surface, which corresponds to a larger ad-atom/surface equilibrium 
distance. The presence of this second stable deep site, even though the energy is higher, 
suggests that the choice of the actual adsorption site is a tradeoff between the chemical 
energy gained due to the larger number of ad-atom-Rh bonds in the deep site, and the 
elastic energy lost due the distortion of the surface when this site is occupied. We shall 
comment more on this issue in the next section. 
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B. Nickel 



In agreement with experimental data [11], for O/Ni(001) we did not find any reconstruc- 
tion at all. In this case the only effect that has been observed is a buckling in the second 
metal layer, whose atoms are not all coplanar but shifted up or down according to whether 
the fourfold site lying above is empty or filled |IT,2"7|-50"|. This behavior is consistent with 
our interpretation of the O/Rh(001) reconstruction: the nickel lattice parameter is quite 
smaller than that of rhodium, and therefore the adsorption site is already small enough 
for the oxygen ad-atoms to bind with neighboring metal atoms without inducing any re- 
construction. For the carbon- and nitrogen-covered surface we have found a reconstruction 
where — contrary to O:Rh(001) and in agreement with experiments findings |||§ — it is the 
filled squares which rotate, i.e. a black reconstruction (Fig. |l]a). In the case of carbon, 
we find that no energy barrier exists between the unreconstructed and the reconstructed 
structures, whereas such a barrier has been detected for N:Ni(001) — though its value has 
not been determined. 



Going back to the results summarized in Table |T|, an inspection of the results for c?oi (the 
distance between the adsorbate atoms and the first metal layer) is particularly instructive. 
In those cases where a black reconstruction occurs (C:Ni and N:Ni), d 0l is small (» 0.1 A), 
thus indicating that the adsorbates penetrate the first metal layer. The penetration of the 
adsorbate into the metal determines an outward (with respect to the center of the square) 
force acting on the metal atoms at the corners of the black squares. This can be alternatively 
be described as an inward force as seen from the center of the white squares. In the case of a 
white reconstruction (0:Rh), instead, the ad-atom stays out of the surface (<ioi ~ 1 A), and 
they pull directly the neighboring atoms towards the adsorption site (rather than pushing 
them away from it) in order to strengthen the O-metal bonds. In both cases, the square which 
is 'squeezed' by the adsorption process (black or white, depending on the case), displays an 
instability towards a rhomboid distortion which is partially opposed by the elastic reaction 
of the substrate. 
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It is worth noting that in the cases where a black reconstruction occurs — and, hence, efoi is 
small — the unreconstructed structure is a metastable (N/Ni) or unstable (C/Ni) equilibrium 
structure characterized by a rather large value of doi (equal to 0.48 A for both N/Ni and 
C/Ni). In order to understand the chemistry of the bond between carbon and nickel we 
have analyzed the local density of states for the unreconstructed and reconstructed surfaces, 
which we display in Fig. || In the two upper panels we show the LDOS projected onto 
the C p x or p y orbitals which are degenerate by symmetry (left) and onto the C p z orbital 
(right). In both bonding and an anti-bonding structure can be clearly identified. 

In the case of the C p x>y orbitals, the bonding structure is resonant with similar structures 
observed in the LDOS projected onto the Ni d xy and d x 2_ y 2 orbitals, while the C p z bonding 
structure is resonant with Ni d xz (middle panels). Inspection of the difference between the 
C p x ,y- and C p z -projected LDOS before and after reconstruction shows that the gain in 
binding energy upon reconstruction mainly comes from a shift of the C p z LDOS towards 
lower energies, thus indicating the establishment of a covalent-like bond between the C p z 
and Ni d xz orbitals. 

IV. DISCUSSION AND CONCLUSIONS 

The behavior of the two surfaces studied here (Ni(001) and Rh(001)) is apparently very 
similar — they have the same STM images and LEED patterns — but a more accurate inspec- 
tion of their structural and electronic properties reveals many interesting differences. The 
most important one is that, in spite of the fact that the adsorption site of the adatoms is 
the same, the kind of reconstruction induced by them is different for O/Rh(001) and for 
C/Ni(001) or N/Ni(001). In the first case the oxygen site is deformed into a rhombus (white 
reconstruction), while in the second case the adsorption sites remain square, and it is the 
empty sites which are deformed into rhombi (black reconstruction, see Fig. |IJa) . Secondly, 
on rhodium the distance of the oxygen atoms from the surface is not appreciably different 
for the unreconstructed and the reconstructed surfaces, while on nickel we have found that 
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FIG. 5. Upper panels: local density of states of the C/Ni(001) surface, projected onto C p x or 
p y (left) and C p z (right) atomic orbitals. Middle panels: same as above, but projected upon Ni 
d xy and d x 2_ y 2 orbitals (left) and Ni d xz and d yz orbitals. Lower panels: difference between the 
LDOS projected upon C p orbitals before and after the clock reconstruction takes place. 
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when the reconstruction takes place carbon and nitrogen atoms become almost coplanar 
with the substrate ones, thus becoming essentially five-fold coordinated. The carbon atoms 
penetrate into the nickel substrate even at zero temperature, and they arrange themselves so 
as to be almost coplanar with the first nickel layer; the nitrogen atoms behave similarly, but 
a barrier is found to exist which gives rise to a meta-stable equilibrium structure in which 
the ad-atoms do not penetrate and the surface does not reconstruct. These facts indicate 
that the mechanism of the reconstruction is related with the atomic penetration into the first 
surface layer: the sites which accommodate the ad-atoms tend to enlarge, and this results 
necessarily in a rhomboid distortion of the empty sites. This distortion has a cost in term of 
elastic energy. We argue that the surface reconstructs when the chemical effects are larger 
than the elastic effects. Upon C or N adsorption, the Ni(OOl) surface reconstructs, while 
Rh(OOl) does not. The reason is probably the different stiffness of the two metals: rhodium 
has larger elastic constants and thus it does not reconstruct. In the case of the oxygen the 
mechanism of the reconstruction is completely different, being due to the O-Rh re-bonding, 
and results in a different reconstruction pattern. 

Further details on this work can be found in the PhD thesis of one of us (DA) 
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